Executive  Summary 


Expendable  Current  Profilers  (XCPs)  do  not  transmit  a 
pressure  signal  from  the  instrument  to  the  ship.  The 
depth  is  determined  either  by  an  empirical  equation  or 
by  comparison  with  a  reference  oceanic  variable.  This 
report  provides  documentation  of  a  computer  program 
for  the  determination  of  a  least  squares  best  fitted 
polynomial  for  determining  the  depth  as  a  function  of 
time,  using  temperature  as  measured  from  a  CTD  at  the 
time  of  the  profiler  drop. 
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I.  Introduction 

A  serious  problem  encountered  when  analyzing  XCP  (Ex¬ 
pendable  Current  Profiler)  data  is  the  error  in  true 
depth  of  the  instrument.  The  XCP  is  free  falling  and 
does  not  transmit  a  pressure  signal  to  the  ship.  Be¬ 
cause  there  is  no  direct  measurement  of  pressure,  the 
pressure  is  estimated  through  the  use  of  an  empirical 
equation  relating  time  (starting  when  the  profiler  is 
dropped)  to  pressure.  Variations  in  launching  proce¬ 
dures,  initial  conditions  and  oceanic  conditions  may 
result  in  significant  departures  of  the  true  pressure 
from  that  estimated  by  the  equation. 

The  XCP  does  return  a  temperature  signal  to  the  ship. 
If  the  true  temperature  profile  (versus  pressure)  is 
known,  then  it  should  be  possible  to  obtain  an  accu¬ 
rate  estimate  of  the  position  of  the  instrument  as  a 
function  of  time  through  appropriate  comparison  of  the 
tru.:  profile  of  temperature  and  that  returned  by  the 
XCP. 

A  program  to  accomplish  this  has  been  written  and  is 
documented  in  this  report.  The  theory  and  some  partic¬ 
ular  results  are  given  In  the  journal  (Oceans,  Septem¬ 
ber  1981)  reprint  in  the  following  section;  the 
detailed  program  documentation  follows. 
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ABSTRACT 

Recent  intercomparisons  of  depths  of 
expendable  bathythermographs  { X8T )  and 
electromagnetic  current  profilers  with  more 
accurate  shipboard  profilers  have  shown  that 
the  empirical  formulae  currently  used  are  in 
error.  A  theoretical  model  for  the  bulk 
dynamics  of  XBl's  (T-7)  and  XCP's  is  proposed 
and  assessed  with  recent  intercomparisons 
between  XCP's  with  accurate  Conductivity- 
Temperature-Depth  profilers  (CTO's).  The 
records  for  the  intercomparisons  were  obtained 
in  the  Gulf  of  Mexico,  where  a  rich  temperature 
fine-structure  provides  numerous  features  in 
the  profiles  that  quantitatively  improve  the 
XBT/CTD  and  XCP/CTD  comparisons.  The  theory 
yields  fall  rate  and  total  depth  versus  time 
predictions  that  agree  with  observations 
somewhat  better  than  the  manufacturer’s 
empirical  formulae.  The  theory  also  furnishes 
a  rational  basis  for  establishing  deployment 
techniques  and  probe  modifications  that  can 
improve  depth  accuracy  and  fall  rate. 


1.  INTRODUCTION 

During  the  past  decade  the  production 
technology  for  expendable  oceanographic 
profiling  instruments  has  improved 
significantly,  consequently,  expendables  are 
becoming  increasingly  important  in 
oceanographic  survey  work.  As  the  measurement 
consistency  of  these  devices  has  improved,  the 
users  have  increased  their  demands  of  quality 
of  the  data  obtained  by  the  probes.  Probe  fall 
rate  and  depth  as  a  function  of  time  have 
received  relatively  little  attention  in 
technical  literature  to  date.  The  expendable 
probes  such  as  XBT’s  (bathythermographs),  XCP’s 
(Current  plus  BT),  XSV’s  (Sound  speed  +  BT)  are 
not  currently  equipped  with  pressure  sensors; 
thus,  the  probe  depth  must  be  estimated  from 
empirical  formulas  based  on  observed  fall  rates 
and  depths  as  a  function  of  time.  These 
empirical  formulas  do  not  in  general  provide 
the  best  possible  estimates  of  probe  fall  rate 
or  depth.  Motivated  by  our  own  needs  to 
improve  probe  fall  rate  and  depth  estimates,  we 
have  ui^e”"taken  development  of  experimental 
methods  and  theories  that  will  help  reduce  fall 


rate  and  depth  errors.  In  this  paper  we 
present  results  of  a  simple  dynamical  model  for 
streamlined  bodies,  such  as  expendable 
profilers,  falling  through  water.  The  model  is 
applied  specifically  to  XCf'  in  this  paper, 
but  we  note  that  it  is  also  applicable  to  other 
types  of  expendables  (Green,  1981).  Next  we 
describe  calibration  methods  that  reduce  depth 
error  by  analytical  treatments  of  XCP  data  and 
profiles  from  reference  instruments,  such  as 
conductivity-temperature-pressure  prof i lers 
(CTD’s).  We  show  that  the  dynamical  model  is 
in  general  an  excellent  predictor  of  probe 
depth  as  a  function  of  time.  We  conclude  our 
discussion  with  some  recommendations  that  could 
improve  XCP  probe  performance. 

2.  THE  DYNAMICAL  MODEL 

A  body  falling  vertically  through  water  is 
subject  to  the  bulk  forces  of  hydrodynamic  drag 
and  buoyancy,  which  ere  in  balance  described  by 
the  following  equation: 

ildv/dt  *£/»w  CDSv2  -  g(M  -  Mw)  (2.1) 

The  bulk  drag  force  is  proportional  to  the 
square  of  the  fall  velocity  (v),  and  the 
buoyancy  is  determined  by  the  difference 
between  the  body  mass  and  the  mass  of  the  water 
displaced.  (A  list  of  definitions  of  the 
symbols  are  in  Appendix  A.)  The  expendables 
currently  in  use  operate  over  a  range  of 
Reynolds  number  (R)  between  10*  and  1(F, 
where  R  is  defined  for  a  streamlined  body  as 
the  product  of  probe  length  and  fall  speed 
divided  by  the  kinematic  viscosity  of  the 
water.  The  kinematic  viscosity  increases  with 
depth  in  the  ocean  and  the  fall  speed  is 
relatively  constant  after  initial  transient 
adjustment,  consequently,  the  R  tends  to 
decrease  with  increasing  depth.  If  the  probes 
were  perfect  streamlined  bodies  moving  through 
quiescent  sea  water,  then  we  should  expect  that 
the  drag  coefficient  (Cn)  would  vary 
significantly  with  depth  and  alter  the  drag 
force.  This  tendency  appears  to  be  attenuated 
by  probe  rotation  caused  by  the  axial 
stabilizing  fins,  which  disturb  ("turbulize") 
the  boundary  layer  flow  over  the  probe  and 
increase  the  wake  turbulence.  The  net  effect 
is  that  the  flow  regime  around  the  probes  is 
fully  turbulent.  The  Cq  of  the  probe  is  much 
less  affected  by  variations  in  Re  in  fully 
turbulent  flow  (Hoerner,  1965),  but  we  should 
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keep  in  mind  that  small  changes  in  Cq  over 
the  range  of  the  probe  can  make  significant 
changes  in  fall  rate. 

The  buoyancy  of  a  probe  also  changes  due 
to  unspooling  of  the  signal  transmission  wires. 
The  mass  of  wire  loss  is  partially  replaced  by 
water  that  is  entrained  within  the  probe.  To  a 
good  approximation  the  mass  loss  rate  is 
proportional  to  the  depth.  Under  extreme 
operating  conditions  this  simple  assumption  may 
be  in  error  (Green,  1981). 

Green  (1981)  found  an  approximate  solution 
to  (2.1)  for  the  special  case  in  which  the 
initial  probe  speed  is  equal  to  the  nominal 
fall  rate  [v(z»«)  =  Vj].  The  mass  slowly 
decreases  and  the  drag  coefficient  slowly 
increases  linearly  with  depth;  for  the  model 
the  changes  are  assumed  to  be  linear: 

M  =  M0(l  +  Az),  C0  =  C0(l  -  Bz), 

z<0;  A,  B,«  1 

The  approximate  sc'-ution  for  fall  rate  as  a 
function  of  depth  is: 

v*  -  VT  [1  +  (A  ♦  B)z]*  (2.2) 

and 

z»  -  VTt  +  (A  +  B)  VTt2  (2.3) 

These  approximate  solutions  clearly  show  the 
dependence  of  v  and  z  on  the  variable 
parameters.  To  first  order  in  t  the  probe 
depth  is  determined  by  the  nominal  terminal 
speed  [VT  =  (2M0g 1  ^CoS  )t] .  The 
coefficient  of  the  quadratic  term  Is  also 
proportional  to  the  relative  changes  in  drag 
coefficient  and  mass  (A  +  B).  The  empirical 
equations  for  probe  depth  as  a  function  df  time 
are  typically  quadratic,  so  we  can  compare  the 
analytical  results  with  empirical  fits  of 
observations. 

Generally  the  effects  of  initial  transient 
adjustment  of  the  speed  on  the  predicted  depth 
cannot  be  obtained  except  by  numerical 
integration  of  (2.1).  Results  of  a  numerical 
integration  of  with  v ( z  =  0)  =  -15  m-s'l  are 
displayed  in  Fig.  1.  The  variations  of  A  and  B 
are  less  than  the  nominal  values  for  an  XCP 
(Table  4),  but  the  shape  of  the  curve  is  the 
same. 

3.  EXPERIMENTAL  ANALYSIS 

Present  expendables,  such  as  XBT's  and 
XCP's,  measure  temperature  -  not  depth  or 
pressure.  The  probe  depths  can  be  inferred 
from  equations  giving  fall  time  as  a  function 
of  depth  (pressure),  but  a  more  accurate  method 
is  based  on  a  mapping  of  the  temperature 
profile  of  an  expendable  on  a  temperature- 
pressure  record  obtained  with  a  CTD  in  spatial 
and  temporal  proximity  with  the  probe  drop.  In 


XCP  Fall  Rate  vs.  Depth 
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FIGURE  1 

the  experimental  results  presented  here  the  CTO 
was  lowered  to  350  dbar;  at  this  time  the  XCP 
was  launched.  The  XCP  passed  the  CTD  before  it 
reached  750  m  depth.  This  method  provided 
near-simultaneity. 

The  use  of  concurrent  CTD  casts,  would,  if 
required,  make  the  profiler  unsuitable  for 
large  scale  rapid  surveys.  Thus,  it  is 
important  to  determine  whether  a  suitable 
equation  for  pressure  (or  depth)  as  a  function 
of  time  of  drop  can  be  obtained  from  a 
relatively  small  number  of  intercomparisons 
with  concurrent  CTO  casts. 

We  are  certainly  not  the  first  group  to 
attempt  this.  Smart  (private  communication) 
has  performed  a  similar  study  by  comparing  fall 
times  with  pressures,  referencing  both  to 
specific  features  in  the  temperature  profile. 

We  have  chosen  to  use  as  much  of  the 
ttmperature  profile  as  possible  with  the  object 
of  reducing  random  errors. 

The  data  set  used  for  this  study  was 
obtained  during  the  deployment  cruise  of  the 
Acoustically  Tracked  Oceanographic  Mooring 
TATOM  '79)  project  (Saunders,  Green  and  Bergin, 
1980).  Concurrent  CTO  and  XCP  casts  were  made 
in  a  region  of  strong  thermal  and  velocity 
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gradients  in  the  central  Gulf  of  Mexico  in 
December,  1979.  The  thermal  regime  was 
characterized  by  a  very  uniform  mixed  layer 
underlain  by  a  strong  thermocline,  rich  in  fine 
structure  but  belonging  to  a  single  water  mass. 

For  each  cast,  an  equation  of  pressure  as  a 
function  of  the  time  of  drop  was  determined  by 
fitting  the  'observed'  pressures  and  times  to  a 
function  of  the  form. 

p(t.M)  .  £  Q„t"  (3.1) 

••o 

The  fit  is  formally  obtained  by  fitting 
the  'observed'  pressures  P0(t0)  to 
p(tp,N)  where  t0  is  the  observed  time  and 
P0  is  the  observed  pressure  from  the  CTD 
cast.  These  are  related  by  the  equations: 

P0  =  P0[T0(CTD)] 

t0  *  t0[T0(XCP )] 

T0( CTD  )  =  T0(XCP )  *•  T0(offset) 

It  is  easily  understood,  that  this  procedure 
will  work  only  if  T0 (CTD )  and  T0(XCP)  are 
monotonic  functions  of  P0  and  t0, 
respectively.  We  were  fortunate  that  there 
were  large  sections  of  the  temperature  profiles 
in  all  the  casts  analyzed  where  this 
requirement  was  satisfied.  The  mixed  layer 
provided  a  convenient  thermostad  wherein  the 
constant  offsets  between  the  XCP's  and  the  CTD 
could  be  determined.  A  typical  plot  of 
interpolated  pressure  versus  time,  with  both 
quantities  referenced  to  the  same  temperature 
profile  is  presented  in  Figure  2. 


PRESSURE  VS  TIME 


The  results  for  the  first  order  and  second 
order  fits  are  summarized  in  Tables  1  and  2 
respectively.  For  both  orders  of  fit,  it  is 
obvious  that  the  initial  offsets  varied  quite 
widely,  while  the  fall  rates  appeared  to  vary 
in  a  narrow  range. 

This  range  may  be  associated  to  the  95* 
confidence  limits  computed  using  the  Student's 
t-distribution  with  24  degrees  of  freedom. 

These  limits  are  denoted  as  and  M~&t 

and  have  been  computed  for  all  the  fitted 
coefficients.  The  nominal  value  for  the  linear 
term  (using  the  quadratic  formula),  4.544  dbar- 
s"l  is  smaller  than  the  lower  confidence 
limit  in  Table  2,  while  the  constant  and 
quadratic  terms,  3.1  dbar  and  -6.75  x  10*4 
dbar-s"S  respectively,  fall  within  the 
corresponding  confidence  limits  and  may  be 
accepted. 

Some  values  for  the  quadratical ly  fitted 
coefficients,  together  with  the  associated 
parameters  from  the  model  equations  are  given 
in  Table  3.  For  the  original  model,  the 
constant  term  must  be  rejected  at  the  95% 
level,  while  the  linear  and  quadratic  terms 

cannot  be  rejected.  After  some  experimentation 
with  the  model,  it  was  found  that  an  initial 
velocity  of  about  15  m  s“l  produced 
acceptable  coefficients. 

It  is  highly  unlikely  that  the  probes  ever 
entered  the  water  at  such  high  speeds,  however. 
The  probable  cause  for  the  large  variation  in 
the  constant  term  is  the  method  for 
establishing  the  zero  time  point  (Sanford, 
private  communication).  This  variation  would, 
for  small  offsets,  affect  only  the  constant 
term  in  the  pressure  (depth )/time  equation  to  a 
measurable  degree. 

When  the  empirical  fits  were  made,  the 
observed  pressure  was  differentiated  with 
respect  to  time  and  this  derivative  was  plotted 
as  a  function  of  pressure.  A  typical  example 
is  shown  in  Figure  3.  The  large  variation  in 
the  pressure-time  derivative  indicated  to  the 
possibility  that  the  probe  was  precessing  as  it 
fell,  causing  an  effectively  changing  drag 
coefficient.  This  was  incorporated  into  the 
model  as  a  drag  coefficient,  varying 
sinusoidally  with  depth.  The  magnitude  of  the 
variation  of  the  drag  coefficient  turned  out  to 
be  much  too  great  if  it  were  to  account  for  the 
same  order  of  the  drop  rate  as  was  observed. 
This  cursory  test  does  not  lay  to  rest  the 
question  of  gyroscopic  effects  on  probe 
dynamics,  since  it  is  conceivable  that 
gyroscopic  precession  coupled  with  the 
surrounding  flow  could  cause  temperature 
anomalies  and  velocity  errors.  This  is  a 
subject  for  future  work.  We  therefore  advance 
the  tentative  hypothesis  that  the  large 
apparent  variations  in  the  drop  speed  are 
artifacts  of  the  XCP  temperature  trace  and 
arise  from  the  restricted  flow  past  the  XCP 
thermistor. 
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DP/'OT  VS  PRESSURE 


5.  REFERENCES 


4.  RECOMMENDATIONS 

In  view  of  the  variations  observed  in  the 
pressure  time  derivatives  with  depth  and  the 
large  variation  in  the  initial  pressure  offsets 
we  make  the  following  recommendations: 

a.  a  signal  should  be  made  available  to 
inform  the  processor  of  the  initial  time  of 
launch  and  surface  impact; 

b.  the  thermistor  should  be  mounted  in  a 
position  (perhaps  near  the  tail)  where 
ventilation  would  be  maximum.  This  would 
provide  for  better  intercomparison  with 
concurrent  XBT  or  CTO  drops; 

c.  the  initial  probe  mass  (M0)  should  be 
controlled  during  manufacture  to  within  a  1% 
tolerance  so  that  variations  in  fall  speed  due 
to  mass  variations  can  be  reduced; 

d.  gyroscopic  effects  on  probe  dynamics 
ind  measurements  should  be  analyzed,  if 
practicable; 

e.  develop  a  new  body  shape  nearer  to 
ideal  streamlined  form  with  greater  diameter  to 
increase  current  sensitivity  and  better 
thermistor  ventilation; 

f.  a  pressure  signal  should  be  multiplexed 
and  sent  up  the  signal  wire. 
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TABLE  1.  LINEAR  FIT  TO  PRESSURE  VS.  TIME 


XCP  Cast 

Qo 

Ql 

r  P 

4 

5.65 

4.63 

0.40 

8 

9.60 

4.63 

1.43 

9 

7.18 

4.54 

2.27 

11 

11.88 

4.61 

2.23 

12 

15.03 

4.51 

2.97 

13 

6.44 

4.57 

0.87 

14 

12.77 

4.42 

2.80 

15 

4.66 

4.58 

0.59 

16 

38.10 

4.66 

1.14 

17 

4.19 

4.50 

1.22 

18 

4.52 

4.51 

1.25 

20 

5.02 

4.53 

1.62 

21 

-0.60 

4.85 

0.76 

22 

8.99 

4.44 

1.26 

23 

6.73 

4.55 

1.00 

25 

10.57 

4.47 

2.29 

26 

8.04 

4.53 

1.30 

29 

3.84 

4.63 

0.40 

31 

9.79 

4.55 

1.71 

32 

-4.28 

4.58 

1.66 

34 

11.27 

4.69 

2.81 

35 

8.03 

4.64 

1.79 

36 

12.29 

4.63 

2.46 

37 

5.00 

4.57 

0.89 

38 

7.99 

4.61 

1.55 

iL 

8.504 

4.576 

<r 

7.464 

0.0883 

T 

3.145 

0.0372 

Mt  0S 

11.649 

4.613 

fix--* 

5.359 

4.539 
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TABLE  2.  QUADRATIC  FIT  TO  XCP  PRESSURE  VS.  TIME 


TABLE  4.  PHYSICAL  CHARACTERISTICS  OF  XBT  (T-7 ) 
AND  XCP 


XCP  Cast 

Qo 

Ql 

Q2 

4 

5.65 

4.63 

-0.141 

xlO-5 

0.40 

8 

7.29 

4.71 

-0.653 

X 

10‘3 

1.40 

9 

1.55 

4.70 

-0.891 

X 

10*3 

1.93 

11 

11.26 

4.63 

-0.128 

X 

10-3 

2.22 

12 

3.85 

4.87 

-0.235 

X 

10-2 

2.02 

13 

3.93 

4.65 

-0.526 

X 

io-3 

0.70 

14 

3.16 

4.69 

-0.150 

X 

10-2 

1.79 

15 

2.19 

4.68 

-0.110 

X 

10*2 

0.57 

16 

32.46 

4.85 

-0.124 

X 

10-2 

0.62 

17 

-0.22 

4.62 

-0.714 

X 

10-3 

1.24 

18 

-1.18 

4.70 

-0.134 

X 

10-2 

1.09 

20 

-1.00 

4.69 

-0.837 

X 

10-3 

1.20 

21 

-0.94 

4.87 

-0.228 

X 

10-3 

0.76 

22 

2.85 

4.60 

-0.849 

X 

10*3 

0.51 

23 

1.67 

4.70 

-0.883 

X 

10*3 

0.63 

25 

6.30 

4.58 

-0.649 

X 

10*3 

2.17 

26 

6.20 

4.52 

-0.313 

X 

10-3 

1.26 

29 

6.38 

4.68 

+0.124 

X 

10*2 

0.38 

31 

4.95 

4.68 

-0.758 

X 

10-3 

1.36 

32 

-12.44 

4.85 

-0. 194 

X 

10-2 

1.47 

34 

1.21 

4.96 

-0.148 

X 

10*2 

1.68 

35 

4.94 

4.74 

-0.638 

X 

ID*3 

1.71 

36 

3.55 

4.88 

-0.144 

X 

10-2 

1.56 

37 

3.28 

4.61 

-0.399 

X 

ID'3 

0.85 

38 

2.03 

4.79 

-0.117 

X 

10-3 

1.14 

XBT (T-7) 

XCP 

Initial 

Probe  Mass 

0.740  kg 

1.480  kg 

Water 

Displacement 

0.161  kg 

0.650  kg 

Initial  Bulk 

Drag  Coefficient 

0.125 

(Hoerner, 

0.35 

1965:  torpedo) 

Relative  Mass 

Change 

(nr+xlO"5) 

14.1 

7.17 

Relative  Change 
of  Cd  . 
(m-'-xlO"4) 

0.6 

1.13 

Nominal  Fall 

Speed  (m  s'l) 

6.5 

4.7 

Probe  Length 

0.215  m 

0.  m 

Probe  Cross 
Section  Area 
(m2xl0'3) 

2.11 

u  • 

t*  3.956  4.708  -8.315  x  10'4 

^  7.323  0.115  7.042  x  10'4 

3.085  0.0485  2.967  x  10"4 

(*+■*•  7.041  4.76  -5.35  x  10‘4 

0.871  4.66  -1.13  x  10-3 


TABLE  3 


M10-5) 

B{10-S) 

to  v(z-0 ) 

Oo 

Ql 

02(10*3) 

7.06 

6.00 

0.35  L5 

3.6625 

4.7007 

-0.979 

10.00 

3.6772 

4.6994 

-1.169 

20.00 

3.7283 

4.6943 

•1.613 

30.00 

3.7963 

4.6877 

-2.017 

6.00 

0 

-2.353 

4.7034 

-0.979 

5 

0.622 

4.7022 

-0.979 

10 

2.387 

4.7014 

-0.979 

0.05 

6.085 

12.450 

-6.863 

0.10 

5.269 

8.797 

-3.625 

0.15 

6.966 

7.181 

-2.282 

0.20 

6.556 

6.218 

-1.711 

0.25 

6.207 

5.562 

-1.370 

0.30 

3.912 

5.077 

-1.161 

T0 

Tq(CTD) 

T(offset) 

T0(XCP) 

t 


APPENDIX  A 

LIST  OF  SYMBOLS 

(Mn  -  M)/(M0  700  m),  the 
relative  change  in  probe  mass 
[Cp  -  Cd(700  m)]/CD  700  m 
Initial  drag  coefficient  of  the 
probe  after  impact 
The  probe  drag  coefficient 
Gravity  acceleration  9.8  m  s"1 

9(1  -  VM0  ) 

Probe  range 
Probe  mass 

Probe  mass  at  impact 
Mass  of  displaced  water 
Order  of  the  least  squares 
polynomial  fit 
Pressure 

Observed  CTD  pressure 

Cross  section  area  of  probe  normal 

to  axis  of  figure 

nth  coefficient  of  least 

squares  polynomial 

Temperature 

Observed  temperature 

Observed  CTD  temperatures 

Computed  temperature  offset 

Observed  XCP  temperature 

Time 

Vertical  speed 
(2Mpg7y^c0S)fc,  the  nominal 
fall  speed 
Depth,  z  0 
Mean  value 
Standard  deviation 


l 

j 


i 
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III.  Summary  Program  Documentation 


PROGRAM: 

PURPOSE: 

MACHINE: 

LANGUAGE: 

AUTHOR: 

FILE  LOCATIONS: 

INPUT: 


XCP-FALL-RATE 

To  compute  emprical  rates  of  XCP's  based  on  simultaneous  XCP  and 
CTD  temperature  data 

UNIVAC  1108 

FORTRAN  V 

Kim  David  Saunders 

Absolute,  Relocatable  and  symbolic 
Elements  -  COOE331*WORKFILE.MAIN/XCP 

Unit  8  -  input  FEB  File  -  XCP  data 

Unit  9  -  input  FEB  File  -  CTD  data 


Unit  5  - 

line  1: 

ISEGX,  ISEGY 

line  2: 

PLIMX  (1),  PLIMX  (2) 

line  3: 

PLIMX  (3) 

line  4: 

ANS 

line  5: 

IDEGREE 

ISEGX  -  seg.  no.  for  XCP  cast  on  unit  8 

ISEGC  -  seg.  no.  for  CTD  cast  on  unit  9 

PLIMX  (1)  -  upper  pressure  limit  for  computing  temperature 
offset 

PLIMX  (2)  -  lower  pressure  limit  for  computing  temperature 
offset 

PLIMX  (3)  -  starting  pressure  *or  T/P  intercomparison 

ANS  -  yes  (for  XCP  temp,  smoothing)  no  (for  no  smoothing) 

IDEGREE  -  degree  of  fitting  polynomial 
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OUTPUT: 


Unit  6  -  Pressure  limits  (repeat  of  input) 

-  Short  FEB  file  headers 

-  Averages  of  XCP  and  CTD  temperatures 

-  Temperature  offset 

-  NT,  N2,  Tl,  T2 

-  Coefficients  of  best  fit  polynomial 

-  Standard  deviation  of  fit 

Unit  20  -  Listings  of  pressures,  temperatures  and  fall  times 

Unit  25  -  DISSPLA  Plots 


NT  -  Number  of  temperature  points  in  input  for 

interpolation 

N2  -  Number  of  points  to  be  interpolated 

Tl  -  Starting  temperature 

T2  -  Ending  temperature 


SUBROUTINES  CALLED:  INITLZ 
XCPRDR 
CTDRDR 
TMPOFS 
TMPCOR 
DPTCOM 
PTINTP 
TTINTP 
FITTER 
PLOTER 
INTRPL 

<<DISPLA  Package  subrout1nes» 
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ADDITIONAL 

INFORMATION: 


1.  The  XCP  FEB  File  must  have  the  form: 

4  variables:  P,  T,  U,  V 

2.  The  CTD  FEB  Fi  ,  must  have  the  form: 

8  variables,  P,  T,  (other  6  are  not  used) 

3.  The  maximum  length  of  a  segment  is  1000  cycles 

4.  Each  drop  is  stored  in  a  separate  segment,  and  the  data 
should  be  interpolated  to  1  decibar  pressure  surfaces,  and 
low  pass  filtered  with  a  cutoff  of  about  8  decibars.  The 
filtering  and  interpolation  may  be  done  via  ZFILT  and 
ZINTERP.  Both  of  these  are  standard  FEB  File  utilities. 


IV.  Runstream  to  Create  and  Run  the  Program 


TfST  P  UN  :  * CP/CTO  CAL IRRA  TION  PROGFAN  -AD  SAUNOFRS  -  AA  7  33  -N090A331 


COOE33 

I 

7 

* 

H 

5 

6 

7 

8 
9 

10 
I  1 
17 
1  7 
1  A 
IS 


UWORKFllEll I.HAKFR/XCP-CAIIP 
mil  if  N  .FOOF33I  •bORFFILF 
RASG.UP  XCP-PR-?,F///?flOO 


AASG.AX  XCP-PR-? 

RPltKPT  PRINT*/ XCP-PR-? 

RHDG.P  IF S 1  RUN  :  XCP/CTO  CALIBRATION  PROGRAM  -No  SAUNOF RS  -  AA7JI 
GPRT.S  b.NAKEP/XCP-CALlB 
iPRT.S  b.M AP/XCP-CALIB 


•a st. i  7n.F///7nnn 
RAS&.UP  7S.F///70 00 
OF  REE  ?! 

RASG.AX  ?S. 

SIRS  ?B. 

-RF-RS  70. - - - : 


16 

BAS6.A 

X  N. 

17 

RFOR.S 

W.MAIn/XCP-CALIB 

18 

RFOR.S 

W  .  I  Nl  TL  7/XCP-C  Al  IB 

19 

RFOR.S 

V.X  f  PROP/ X CP— CAL  IB 

?n 

RFOR.S 

W  ,C  IBRD  r/xc  p-cal  IB 

- 

RFOR.S 

6.1  MPOF  S  /  XC  P— C  At  ff - 

7? 

RFOR.S 

6.1  HP  CO  P/XC  P— CAL  IB 

73 

RFOR.S 

6 .DPTCON/XC P— C AL IB 

7  A 

RFOR.S 

N.P IINTP/XCP-CAL IB 

?S 

RFOR.S 

b.  T  1  TNT  P/XCP— C  Al  IB 

76 

RFOR.S 

W.F  IT  IE  R/xCP— CAL  IB 

77 

-  RFOR.S 

W.PLOTE  B/xCP-CAE  TB - 

78 

RFOR.S 

KIISFILF  .TNTRPL/AKINA 

79 

RADD.Pl  N.MAP/XCP-CALIB 

3D 

RPACK 

N. 

31 

RPRFP 

b. 

37 

RPRAPT 

PRINTS 

33 

XCP  —PR—  ?  ».P  6 

3  A 

RFfcEF 

XCP-PR-? 

RPR  T*  S  N.MAP/XCP-CALIB 


V.  Map  Runstream 


TEST  PUN  :  XCP/CTD  CALIBRATION  PROGRAM  -AO  S  AUNDFRS  -  XA773  -N0R0A37| 


COD E3 3  UNO  RAF  ICE ( 1  I .NAP/XC P-CA LIB 

I  RMAP.IS  .  It  .X®  T/AC  P— C  Al  IB  - 

7  IN  W.MAIN/XCP-CALIB..INIH 7/ XC P-C Al IB . . XCPRPR/X CP-CALI B 

-  7 - IN  W.eTDRORFXCP-CAtfB - 

A  IN  N.TNPOFS/XCP-CALIR..TMPCOK/XCP-CALTB..DPTCOM/XCP-CALIB 

5  IN  M.PTIMTP/XC  P-C A l IB  .  .  FT  I  NfP/XCP— CRL  7B  nPITTFR/XCP-CAl IB 

6  IN  N.PLOIrR/XC  P-C Al IR .C  PDF  13  URDSF ILE.7REA0 

7  IM  C8BFM»*«BSFIIE  .TNTRPE/»*t** - 

A  LIB  OISSPl AALIBB-D. 


RPR  I .  S  ACP  -CAL  IP 


-NOROA331 
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VI.  Listing  of  a  Typical  Program  Runstream  with  an  Example  of  the  Output  to  Unit  6 


IfiT  7JL  s  *Cf/ClC  C»lIPr*tIOH  PROOF**  -  *  3  TAuNOERS  -  X07I3 
COut  JM.-oPRflLe*  1  ».RUN-l/MCR-C»Lla 


1 

•EPEE  25 

y 

AASG.UP  FCP-2S.E/ //2000 

5 

•Free  XCP-2S 

4 

-USE  25, XCP-25 

s 

o*SS,AX  25 

6 

•»SG,TVJ  IN1,U9Y,S937 

T 

«ASG,T  9.E///2000 

4 

<a»SG,T  8,F///?000 

9 

•POWE  IM.U 

1) 

WCOPT.G  INI, 9 

11 

•Eree.s  INI 

u> 

-*S6,TVJ  I N?  ,1)9  V  ,  0&  39 

11 

OCOPV  IN?  *  8 

10 

-tree  in? 

IS 

«*0T  W.XRT/XCP-CAL1B 

16 

<*.! 

IT 

60,73 

14 

120 

IT 

T  £5 

23 

1 

2 1 

•  ERE  t  25 

■A&G.I  20 

,F///203Q 

facility 

WARNING  1GOOOCOOOOOO 

•ASO.JP  XCP-2S.F///2000 
FACILITY  BEJECTEO  000300033030 


•  f«el  X  CP-25 

FACILITY  warning  lOOQOPOQOQOn 


XCP-25. 

ready 


AERS  XCP-2S. 

F URPUR2 7R3B  £33  RLT3R1  00/10/81  15:18iST 
END  ERS. 


*USc  25,XCr-25 

REACT 


•E«S  ?0. 

F  UMPUR  Z  7R  38  £33  RLT3RI  00/10/81  15j!9j00 
END  f RS • 


i»SG,Ai  y. 

facility  warning  lozaoouooaoo 
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-NOPOA J Jl 


Tf  ST  RJP  !  *Cf>/C?0  C*lIg#*TIQN  PROs«.A*  -  k  Q  SAUNO£»S  -  XAT33  -KOOO*iIl 

•ASG  TX J  In1,u9»,S9J7 
9  EaO» 


V,f///?000 

fdcuifr  karnins  lruooojTDoan 


•  ASG.T  8.F///7000 

facility  mtkiNS  louoQcooaaoo 


j*G«  IM,0 

FUt.PCi»?7R3»  f33  R  L  73*?  1  09/1C1/P1  ISSIRSOT 


•COPY.S  INI, 9 

CO  bf  J 31  *OU  TF  1  IQ  I  COPIED  ON  06/11/90  *  T  lb:??:S3  SPTUOS  REFL-5R3TK 
;-9  CLOCKS  CJPIEO. 

[Of  EVCOUNTEPrp  ON  INPUT  TAPE 


■FPLE.S  INI 

react 


•ASG.TYJ  I H2 f U9X f 0  639 
Hi  AUT 


«CUPT  In/, 8 

FURPUR/7RJ8  E  33  RL7JRI  09/10/91  1 5 j ?0 t 5? 

839  PLOCKj  COPIED. 

EOF  ENCOUNTERED  ON  In*UT  TAPE 


»F*El  X  *2 
PE  AOT 


• *07  9.*0T/*CP-CAlIB 

enter  the  Input  seepents  for  the  *cp  ano  ct3  data  in  free  fopnaT 


ENTER  Tme  UPPtR  aNO  loner  PRESSURE  lIMtj  For  C0*PuTlN6  THE  TEMPERATURE  OFFSET 


PLlr*(l>  -  S3. oooouo 

PLIPTI/I  i  70.000303 


enter  starting  prcssure  for  t/p  interconf«rison 


R&  F  9X-SP  9  C  R  ATO  ‘  IS  IS  *0 

00  TO J  PANT  SHOOTHINS  OF  «C®  TENREb*TUPE5  ?T 
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1 1  v, |  :  *cp/cjr  calibration  fro^tan  -  ►  o  saunofrs 


X4TJJ  -NoP£)A5Tl 


Bu  p  7U  3c 

T*av  = 

TCAV  I 

t  c«rsrt: 


1  10303® 

?S.41 1 
2*.. 473 
*  •>  1866  ~ul 


<*9j  r  t»o  i*o  ion 


VT  :  577 

*2  X  126 

II  :  24 •  c>0(i03 0 

t2  =  u  .  rroooo 


EMIT  TmE  DtrPEo  or  flTT]*^  POL»NO»1*L 


1  «II> 


1  s.<rl4*4f 

2  4. 374780 

SI  ANOA  PU  OE  w  IAT  I  ON  OF  PMfSSU.H  =  .91*»t>JJl4 


81014 " OF  DISSPLA  fLOT  G£  Nf  RA  T I  ON • 


PLOT  40.  1  .ITU  THE  TITLE 

OP/OT  IS  PRFSSORE* 

NAS  4E  F  4  COIPLETEO. 

PLOT  10.  RE*oS 

plot  i  is. 21. as  fri  to  apr,  19« 1  joh=eu*o5  ,  naiocea  oisspla  her  »-n 


0*7*  roo  plot 

40.  or  curves  o*amn  1 


N0RJ7.  AVIS  LENGTH  7.J  INS. 
»LPT.  tits  LtNBTH  9.0  INS. 


0417.  ORIGIN  .3013  »t»T.  ORIGIN  .0300 


MORI 7 *  AUS  LINEAR 

STEF  S 12C  .I429»ci  UNITS/INCH 


13 


ffiT  pu*.  s  *cp/cTf!  calinRaTiOn  reon***  -  »  o  saurocrs  -  xrtjj  -nopOaui 

,tCf.  **1S  tlNtAR 

Sirr  SIZE  .bbb  7»02  UNiTS/lNCH 


LOCATION  OF  CURRENT  PHYSICAL  OR  161*. 
*:  I. 00  f:  .SS  INCHE* 

from  loner  left  Corner  of  pa6£ 


test  RJN  :  XCp/ClO  CALIBRATION  PP0SFAM  -AD  SAUNDERS  -  XA7jJ  -N0R0AJ5I 


°l.C  TT  I  n6  COMM£*ClNf- 


.  DISSPlA  VERSION  8.0 

NO.  OF  FIRST  PLOT  <£ 


PLOT  NO.  .  mITH  TMF  TITLf 

OP/OT  VS  PR r SSORE  -  FILTERED* 

M*S  3EE4  COHPLETEO. 

PLOT  ?P*  "Vs^O.ON  FBI  10  APR,  19R1  UOOrEUXflS  t  '"AROCEA  OISSPLA  »ER  8.0 


DATA  FOR  PLOT 

N} .  OF  CURVES  OR  ANN  l 

HORI 7.  AVIS  LENbTH  T.O  INS. 

VERT.  AXIS  lEN&T H  R.U  INS. 

MORI? •  OR  16  IN  .3113  VERT.  ORIGIN  .0300 

MORI?.  AXIS  LINEAR 

STEP  S I?E  .IN2R*01  UN|TS?INCH 

vert,  axis  linear 

STEP  S 12c  .bbb T*0?  UNiTS/IMCM 


.  LOCATION  OF  CURRENT  PHYSICAL  0RI6IN  . 
.  nz  1.30  *=  .SS  INCHES  . 
.  FROM  LONER  LEFT  CORNER  OF  PAGE 


CNO  OISSPLA  —  29  Jo  VECTORS  SE  NCR*  TEr  IN  2  »LOT  FRAMES. 


IGNORED 


-  IN  CONTROL  NODF 
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*  NJ 


Program  Listings 


Tf  S  T  PUN  :  XCP/CTO  fALIRRATJON  PRCROAN  -  K  f)  SAUNDERS  -  XATIi  -NOPDA33I 


COOE33 IRVOPNt ILF (I  I ."A IN/XCP-C  ALIfi 

I  £*tt  HUH? 

CAU  XCPRDR 

- tAtt  CTORCR 

A  CALL  THPOF  S 

5  tAtt  HtRF» 

6  CAU  OP  I  COM 

7  CAtt  OTTArlP 

R  CA U  IT  I RlP 

- eAtt  niio  - 


ir 
1 1 
i? 


CALL  PlOTFR 

-  C-1AU 

J  fwr 

t  NO 


■  }MHIAtL?f  PRB6RAH 
i  RF  AD  XCP  DATA  FPOH  FFR  Flit 

*  CALCIIL  A  TF  TCMPFRAIURi  OF  F  St  T 
O-  COAMAf-R-t  -TACO  TFA*Pf  »A  T  MRf  S 

*  SF  L  r  CT  COMMON  OFPTH  INTERVAL 

*  H»ffPPDt*H  PRESSURE  ON  TfMPERATURF 
d  TNTF  FPOL  A  If  IINF  ON  TFNPFRATURf 

a  PLOT  THE  PRf S  SURF /TIME  CURVF 


1 


4 


lest  RUN 


XCP/CTO  CALIBRATION  PROGRAM  -  K  n  SAUNDFRS  -  XA733  -NOROA331 


COOE331*NO»KFUF  (1 
I 

?  c 

—  S - € - 


I.INIT17/XCP-CALIP 

SttMRWlT  I  Nf  t  N  t  tt  7 


A 

5 

6 
7 
A 
tt 

10 
11 
l? 
I  t 
1A 

15 

16 
17 
16 
19 
?o 
f-h 
7  7 
7  3 
?A 
?5 
7b 
7  7 
?A 
79 

30 

31 

32 

33 
2A 

35 

36 

37 
3A 
39 
AO 
A  1 
A? 
A3 
AA 
AS 
A6 


C - 

C 

C  PUR 
C 

C - 


POSf  :  TO  INITIAL T7F  THF  XCP  CALIBRATION  PROGRAM 


700 

20U0 


COMMON  /  01  AGS  /  lOIAGttO! 

COMMON  / ROATA/  WRI 100001 
COMMON  /ROOC  ]  /  TBOCRIlOOl 
COMMON  / RDOCF /  FDOCRIlOOl 

COMMON  /RHOR/  IRHORIinnt 

common  /blocki/  tiMfxttootii  .txtiooftt.pxt  tonoi.tct loom. pci  loooi 
common  /block?/  plinxuoi.piimci  ioi  ,ise6c.i SE6x .tofsft 
COMMON  /BLOCK  3/  TEMP! 10001 .P 7*000  7, TIME I lOOfll 
FCAl  BUFFER  110. lOnOI 

oRITMG.  1001 

Hr«ois.?uni  jsfgx.isfgc 

FORMATI*  ENTER  THE  INPUT  SEGMENTS  FOR  THE  *  . 

•*CP  AND  CTO  BATA  IN  ROE f  FORMAT •//! 

FORMATI 1 

- NRI  IF  1 6  .401* - - - 

REA  01 5 ,  700  I  PLlMXll  1.PLIMXI7I 
hRITflb.HlTI  PLlMXll  F.PL  TMKI7  7 

FORMATI*  PLlMXll I  r  *,G7n.*/*  PLIMXI?!  r  «.&?0.t//l 
IF  1  PLIMXIH.Lt  ,P1  (MX  IP  IE  GOTO  IftOO 

IF  1  PL TNXtl I  ,F0.  PLIMXI7I  1  GOTO  7000 

—  — P  '  PL  I  MX  H  7 - - - - 

PL  IMX1  1  I-PLIMX  171 
PL  INK  (7  I-  P 
CONTI  NUF 

FI  INCH  I  =  PLIMXllI 
PIIMC171  :  PLIMXI7I 

FORMAT!  *  ENTER  THf  UPAfS  AN&  LOUFR  P»F CLUCK  LIMITS  For  * . 
'COMPUTING  THF  TFMPFRATURF  OFFSET'//! 

WR|Tf  16.  10?  1 

HF  AO  15.200 1  PIINXI3I 

PL IMCIII=PLIMX|3| 

I ORM  All*  FhTFR  STARTING  PRESSURE  FOR  T/P  INTF  RCOMPARI SON  */ /! 

RF  TURN 

LND 


APRI.S  W.XCPROR/XCP-CAlTB 


TEST  BUN  :  XCP/cTp  CALIBRATION  PROGRAM  -An  SAUnDTRS  -  XM7Y5  -NOBftATTI 


coot  33  IaWORaF  I 
I 


Lf»l  I.XCPRDR/XCP-CAL IP 

bUBROOT  iNf  ACPPftR 


- .  I 

PURPOSE  :  TO  REAR  IN  XCP  DATA  FROM  LOGICAL  UNIT  b  INTO  THE  I 
M*  I*  ftUFf  FR  ***** — ANft  TRANSFER  THE  OAT* T»  TW  I 

PRESSURE  ANP  TFMPRATURE  ARRAYS.  TO  COMPUTE  THE  XCP  I 

- TE  HE  ARRAY  FROM  fitf  PRESSURE  ARRAY.  USING  I  Mg  - 1 

INVERSE  OF  Ton  SANFORD'S  ALGORITHM.  T 


- COMMON-  APIA6S/  WAGHW - 

COMMON  /ROATA/  VRflnOPOl 

r  a  mm  AAi  j  poor  t  i  I  ftftfa  i»nnj 

vOTnwW  F  "1IW  v.  I  T  1  Uw  Vlrtl  U1T  f 

common  /poocf /  Foocpimni 

COMMON  /B0OCA/  «OOC«4  TOOT- 
COMMON  /RHDR/  IRHORIIOOI 
COMMON-  /Bt  OCR  !-/-  TIME  X  < 


i  .PC I  |«WO* 


COMMON  /BLOCK?/  PLTMX  HOI  .PI  IMCI  10|  .ISFGC.I  SEGX  .TOE  SET 

common  /block**  tempt  tnnn i  ,p « i nno*. ttmei  Term* -  - 

COMMON  /BLOCK  7/  VES.AMS 
PFAt;  B0FFFR  C*  .  1  ft  Oft  4  »*.ft  .f.T.Xl  IOBOT  — 

INTEGER  yes.abs 

f  OtHEEOfc  f  NCF  4  YRtTnAHf  F  f  » 1 4  .  f  ♦  ».  «  VO  f  WTO  HtHHW - 

YFS  r  nhYES 

IDIAGI?)  X  J 

IDIAGI  it  =  1  ft  Pfl  - 

ioiagiai  x  iunn 


I01AGIS1  - 
IDIAGIbl  X 
101*64  44-s 
IDIAGI A  1  x 
1ft  |*ft<  At  X 
IftlAGI  in  lx 


+0 — 
MO 
*0- 
mr 
> 

1 


0*4*  F.0»**4.-tA*.S*A. -0.000b  TUBE 


f  i  c  *n  tur  R  I*rr  f  P 


4*0 


UO  AAA  I 
VRIII  x  ft 
C6MT  TNUf 

JO  9  AH  I  x  J  .  1000 

- 0*04-  A- ft - 

T  XII  I  x  ft 
T-fOf  X  4*4x0 
CANT  1 N II E 
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TfSI  PUK  s  XCP/CTD  CALIBRATION  PROGRAM  -mo  SAIINOFRS  -  XA7T3  -MOROA331 


C - 

C - 


no  loiio  i  -  t.tono 

IP*  -  WfPW  1.11 - 

i ipx  .6t.  i nnn  .or.  ipx  .lt.  i  i  &oto  mno 

- WFHWH  =  RUFFrRtt.lt - 

TXlIPXI  =  RUFFFR 17.11 

TI*CX<  IPX  m-mWFAAF^Utf-tfmWFjAAt 
CONTINUE 


C 

C  SAM 
C 

C 


TN  XCP  TFMPFfiATURFS  IF  DFSIRFD 


•RITF 16.1001 

«F  AD  18.1  011  AD8  - - 

FORMAT! •  00  TOO  WANT  3M00THIN6  OF  XCP  TEMPERATURES  ??•//! 
f ORMAI 4 1  *6  1  -  _ _ _ _ _ _  _ _ 

IFIANS  .  Nt .  YfSl  RETURN 

M  HiflF-I  -  3.998 - - - - - 

1 X  XI  I  )  =  .75*  ITX  IIM  l*TX  I  I~1 1  1*0  .?*<  TX  1 1  «?I*TX  1 1 -7)  l«n.3«T  XI  1 1 
(MIIME 

00  1007  I  =  3.99  8 
3X41  l-TXXIll/}.? 

CONT  1NUE 


APR  T. S  W.CIDRDR/XCP-CAilb 


I 


TfST  RUN  :  YCP/CTO  CALIBRATION  PROGRAM  -HP  SAUNPFRS  -  YNTYT  -NOROA33I 
COOF3 31*kORKF iLfl 1  I.ClDROR/XCP -CAUP 

I.  t  iin«\  Am  X  I  AiC  X  r%r%  t%  n 

juB"  ini  T  J  Ht  L  T  IM u "  —  - - 


TO  RFAO  IN  f TO  OATA  FROM  LOGICAL  UNIT  8  INTO  T MT 
OOfFfP  ARRAY  AMP  IRANSftR  YHf  OAT  A  TO  T+4f 
PRFSSIIRF  ANT)  TfMPRATORf  ARRAYS  •  TO  COMPUTF  THF  *CP 


INVFRSf  Or  TOM  SANFoRO’S  ALGORITHM. 


C - - - 


COMMON  /ROATA/  YRTIOOOOI 
COMMON  /ROOCI/  fOOCOUOOT 
COMMON  /RPOCF /  FOOCRIinni 
COMMON  /VOOCAC  AOOCMCHTOT 
lOHMON  /RHOR/  IRHORT ICO  I 

common  /plock?/  pl  i  my  <  ini  .pi  imc  i  tni  .isfgc.i  sfgy  .tof  sft 

COMMON  COLO  CM  3/  TFMIH  1THHH  .PTYOOOL,  THtFL  LOflOT 
COMMON  /PLOCK  7/ VC S . ANS 

in  if  of p  YFS.AOS 
FF  Al  HUfFFR  C8.10DPT  .TCY  urnni 

f  OUIVALT  NCFlYRd  I.BtlFim  I.lll 
F-OLMVOLFNCI  <  YB  ♦  H  .MTFffBLT. in.  1¥» 41001  I.HKIU 

lOIAcm  r  1 

1PIAM7I  =  I 

iniAMii  r  in  on 


lPlAulbl  =  ion 
H*IA6(T»  :  ion 
miAbiBi  -  mr 
lOIAOlN)  -  i 
HIT  AC  I  in  1  r  l 


t  — 
C 

C  CLC  A 

C 

C 


R  THf  BUFFI R 


BO  ROY  I  -  LAtonon 

YRII )  =  n 

CO  NT  |N  UF 

iio  *>*»*  1  =  1.  mnn 
pci  I »  =  n 
T  c  ( 1 1  =  n 

f  AAtf  14|A|r 
V.  * r“  *  l>TUJ 


9*8 - 

t 

C 

t  L  CAP 
C 


IN  IMF  CTO  BATA 


»♦*»♦♦»»»»»»»*»»»»»»♦»»♦»» 
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Trsi  RON  :  XCP/CTO  f 61  1RR4 T  10  N  PROE.FAH  —  44  0  S64INOFRS  -  X9733  -N0RD4  331 


<5  7 
SR 
S9 
bO 
61 
67 
63 
66 
65 
-6b 
67 
6R 

69 

70 

71 
77 
73 
79 

75 

76 

77 
7ft 
79 
RP 
R  1 
R? 
ft  3 


7007 


2010 

loon 


C*ll  7Rf 40|9  .If .ISF6C1 


no  inoo  i  -  l.inoo 

IFHMWttl  .Nf.M  6670  7006 
IPC  =  HUFFFR(1.T<  ♦  .6999 
IFIIPC.LT.1  .06.  IPC  .6  7 .1-660  1  6676 
PCIIPO  =  nuFFEMl.il 


-  Eftift  a non 
V«  I  V  “  T  “t 


iooo 


CONTINUE 

IF (166064 1 1 .NC .161  6076  7010 
JPC  -  BUF10ll.il  «  .9999 
IF  1 1  PC  .1 T  .1  .06.  T  PC  .61.1000!  6070 

pciipci  r  bufi n • i. n 


1000 


fl  5 
R6 


CONTINUE 
CONTINUE 
IF  I* NS  .NF.  TFSI  RETURN 

t,o  inni  i  =  3 .99  r 

TCXC  II  =  .7S*|Tril«l  1«TC1  I 'll  I  ♦0.7941C  4  I*?»«TC4  I“?  1 1*0.  39TC4  II 

7601 - tONT  INUF  -  - - — - - 

no  inn?  I  -  3.99  8 
?C(I»X?Cx4I>/l.7 
CONTINUF 
6FIU6N 

NRITECb.mi 

- FOR*  614  *-  66606-7*  ffH-  FIEF  &76UC  TU6F  —  If  6NIN4T1N6*  I 

STOP 

i  NO 


104,2 

9999 


4PRT.S  b.TNPOFS/XCP-C*Ll» 
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test  r un  •  xcp/cio  calibration  program  -bo  saunders  -  *a7tj  -nordajti 


COOEjn«NO»KFIlf  I  1  I.TMPOES/XCP  -C*l  IT 

I  SOBROtiT  It*€  TMPOF  S 

?  C 

-  e - - - - - - — 


C  PURPOSE  :  THIS  SUBROUTINE  COMPUTES  THE  AVERAGE  OFFSET  T  E  MP£  R  ATURt 
C  OXf  B  A  Biff*  ASStmfU  OfPTH  RANGE  f  BBTAIMEft  H*  UrtllPt 

C  FOR  CORRECTING  THE  XCP  TE  HPERA  TURf  S . 

^ - - KVfE^E  I  ^  AVERA6EL  TCTB— WtCPT - 

t 


COMMON  /  G I A  G  S  /  I  01  AG { 1 0  I 

camm aai  icnit  a  a  uni  i nnnn  > _ _  -  _ _ _ _ _ _ 

■  ttfWWlflr  f  WITH  ?  H  “  TTM  W 

common  /Rnoci/  mocPiinni 

COMMON-  /HOOCF  /  f««»EHW(  - 

COMMON  /R00C4 /  AOOCREIOOE 

r  a  mm  At.  j  mo  n  m  1  nuftn  i  i  no  t 
V  tT T/rT  ^  FfFTtFff  /  *  Ifni/  FT  f  J  ETTT  / 

roMMON  /BLOCK  1  /  TIMEX  (10001  .Til  1000  |.PX  (  1 00  01  .T  C  I  10  001  .PC  (  |  01  0 » 
fOHKPH  /BL-OfB  P/  PL  tMELI-ftt  .PL  Tuft  LEH-tfSEGE  .tsEGX  -«  t«P  VE  T 
COMMON  /ML  0CK  3/  TE MP( 1000 1 .P( 1000 |. TIME  | 103UI 

Gf*L  EHTEEf«I+&.t«H>I  -  - - 

FUUlVALFNCr  (VR(  II.RUrr(R( I .1 1  I 

NfAt  TNAV.ICAV  -  ■  -  - 


TX  AV  :  l) 

LCAR-s  « 
n  ;  n 

It  -  PL  I  MX  I  t  I 
IP  -  PL  I  MX (01 


com mon  /block  J/  Tr*P« inno» .p«inno».ilnM iono> 


rouivALfMCC  (vpi  ii.nuFFiPd.m 

ft  =«rM**<3»  ♦  o^ont - 

I?  -  II 
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TfST  RUM 


XCP/CTO 


CAl  IRRAIIOn  program  -mo  saundfrs  -  xatti  -noroaiii 


COOf 331*NORKFILf II  I 
1 

7 

4 

4 

6 


.PT  JNTP/KCP-CALIP 

MjBROUTlNf  PTINTP 
COMMON  /OIAGS/  IDIAGdOl 
- COMMON  /ROAT4/  4P4100B01 

common  /woci/  rnocPJtom 

COMMON  /ROOCF/  FBOCP4MH>» 
COMMON  /POOCA/  AD0CRI1DUI 


COMMON  /OHOB/  IP NON 4 )OA  A 

common  moon/  timi *  i loom  .ixuooo  i.p* i  innm  .ic  t  loom  ,pc<  mnoi 
PL  f MX  ClfH  .Pf  4*0  101  rfWtt,litH  .INF&fT 


1  n 

COMMON  /FLOCK  3/  If MP 1 1000  I .P 41000 1 . 1 INf 1 

1 1 

COMMON  /PL0CK4/1FMP7<  10004  .07410001  »HMf  ; 

1? 

COMMON  /RLGCKS/T 1.I7.0T .N1 .N? 

l  i 

pFAi  RufFf  r  «m.i  non  i 

14 

[OUlVALFNCf  <  WR  (  |I.RUfPfR<l.lll 

!«. 

Of  *  <f»l  -  -  - — - 

lb 

1 1  -  pcimx  i3i  ♦  n.nm 

1  7 

17  -  PllM*  »4»  ♦  0.001 

IB 

r i  =  him  t 

|4 

17  :  T  X4  1 7  1 

70 

if«Ti  ,c,t.  toiiii  n  -  tcnn 

7V 

41  4-17  .44.  KCI7H  -»7--s  4C+W 

77 

IT  1  =  10*1  1  -  0  .S 

21 

IT?  ;  10* f 7  ♦  0.4 

74 

fl  -  ffl/IO. 

7S 

17  =  117/10. 

7b 

no  iooo  1  ;  11.1? 

7? 

-  4-4  i  4 - H*  4  -  -  — 

7* 

IJ  :  12  -  ix  n 

7<J 

TfMPIHl  ;  1CIUI 

in 

P«If»  :  PCIIJI 

11 

inco 

CONT INUt 

17 

NT  -  17-  11 «  1 

11 

N?  -  -0  .  . . 

14 

no  mm  i  ;  l  .moo 

}*. 

N  7  -  N?  ♦  1 

lb 

TfMP?IN71  r  T?«0T*«I-11 

17 

lf«UMP7(N7l  .Gf.  in  GOTO  100? 

1  A 

inr.i 

CONT  INUf 

14 

10U? 

4.7  =  N?  -  1 

4P 

no  i  util  i  =  l  .Iona 

4  1 

p  7*  1 1  =  n 

*7 

moi 

f ONT INUf 

«  » 

CALL  INffiPUb  .NT  .TFMP.R.N7.TF MP?  ,P?4 

44 

HI  TURN 

11 

f  NO 

RPPT.S  N.T TINTP/XCP-CAl  Ift 
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TEST  RUN 


XCP/CTO  CALIBRATION  PROGRAM  -  A  0  SOUNDERS  -  XATSJ  -NORDAJS1 


COOt  33  1AWORNF  IlE  I  1  1  .  T  T  IN  IP  /XCP  -C  Al  I  f 


I 

7 

V 

A 

5 

6 
7 
P 


SUfMMHJTJNf  1 1  INTI* 

common  /niAfis /  roiAfiimi 
common  /tww-mtmwiH 


COMMON 

COMMON 

COMMON 

COMMON 


/rooci /  inocR lino i 
/  roocf /  fooeoimoi 
/rooca/  AoocRimm 
/«mor/  iohori  mrrt  — 


in 
1 1 
17 
I  I 
1 A 


TIMEX ( 1Q001  .1X11000 I.PXt 100  01 .TCI  10001 .PCI  10101 

PtlNM  tfllrPL  WWWMAWtHWaHHH - 

COMMON  /BLOCK  1/  TEMPI  10001 .PI  1 000 1 . TIMF I  10001 


COMMON 

COMMON 


/BLOCK  1/ 
/BLOCK  t* 


lb 

17 

IP 

l» 

70 


COMMON  /PLOCKS/T 1.T7.DT .NT  »M7 

Rf  At  OKfFFR  l»O.I  «M1  *  - 

rouixAiTNCf  (xeiii.Burrrpu.nl 

- H  -  pl I mk in  »  n^wi - 

17  r  PL  I  MX ( A |  ♦  n.noi 


«o 


lone  i  z 
11=1- 
1  &  =  ■  17- 
TFMPI III 
T1MCI l»> 


»I? 


I  I. 

II  ♦  I 

IT-  A-  1- 

;  IX |I Jl 


-g  T  T 


7? 

?1 

7A 

?S 

76 

77 


inoo 


3  "0 


CONT  INUE 

-CALL  iMTOPLCM-^tT.TfMP^TTMf  ,N7»TF.»P7.T  TMf  7)  - 

ARlTFI6.TnniNT.N7.tl ,T? 

f  0 RH  A T  f  t  //  *  t tt  lift  *  / 

•  NT  -  *.15/*  N7  f  * .IS/*  T 1  =  • .670 • B/ •  T7  =  *.6?0.B//1 
m»|T|I?0»1O0T - - — — - 


7  P 


30 

31 
37 

-T-3- 


inc? 

loo 

2PG 

-e - 


titi  i no?  I  r- 1 .  inno 

^mn?n-.70BiTir»(lKT»ii»>iiw«iiitrciH>Ttm 

CONT  TNUC 

FORMAT*  |H|  .»  PlSPt  AT 
f  ORN  ATISX.IS.LG70.6I 


3A 

35 

36 

37 
3  P 

-30- 


.Rl?r  ipmooc 

no  inoi  i  -  i 


mm 


1000 

ASK  c  370.  |  PPL  fHCMO/l  It »P7I|  »,IFMf?l  O  - - - 

CONTINUr.  *  *  ,  -t  > 

f  0OM  AT  I  H*§  ♦  ■  niSPLAT  OP  ABORTS  Tf  MPp  »Pp»  TIME  >*//-/-> 


a  n 

At 

A7 


iro 


FCRMATISX.IS.7G70.ei 

srwwt  -  — 

r  Ni) 


TfST  RUN  :  XCP/CTO  CALIBRATION  PROGRAM  —AO  SAliNDFRS  -  **7IJ  -N09DATTI 

rOOfS3l»NO«AFtir«l I.FlTTFR/XCP-CAUR 

1  SOfiROUI+Nf  F !  «F«  -  -  - 

7  COMMON  /OI AG  & /  IP7A6I10I 


COMMON  /ROOCI/  IDOCRUnni 

COMMON  /ROOCF  t  FBMAFWF-  - 

COMMON  /ROOCA/  AOOCRdnnl 

COMMON  /BH6M7  IP  HOB  4  1001 -  ' 

common  /plocai/  timf x unnoi  .tx  1 1000  i.px  i  lonoi  ,ic » loom  .pcm mn» 
COMMON  /4>L0GM»A-  PIIMAHOl  .PI  IMP  tOI  .INFOC.I  Mtl  w4QF-*FC  - 
COMMON  /FLOCK*/  TFMPUOOOl.PClOnOl.IIMrOOOQI 
COMMON  /AtflOU/  IT  MR  71  .p?i  mooi.iiMiRtioooi 

COMMON  /BtOCKS/T 1.T7.DT .NT.N? 

Mf  At.  Ntt0004.C0f  F!lfTfV04.U4MTOO»l.I*I?l  IOOOI.I T 1  I!  10001 
RFAL  ALPHA!  1 000  I  ,8M  A  1 1  OoO  I  .M  i  1000  » 

fOUIVALfNCF — I  MB  I  It.  ALPHA!  >  >  >  .1  RR  t  1001*  y*F  TA  i|  UUWtWtt. 

TT II li  II. 

IVR  ITOO  I  I.T  II  ?!  |  It. IVRlNAOl  1.114  Ml  ll.INRI%001l.AI*  II. 

(NR  I  fcOO  II  .C  OFF  I  |  II  .  1XRI  TOO  I  I  .Hi  1  »  I 
f-O  I  GOO  |=  1  •  WAfi 

win  -  l.n 


hr  ITFI  #>.  jnn» 
k(  ACiS.JOl  .EkO-MSi  N 

IF (A  .FO.  n  .OR.  N  .61.  Rl  Rf IUR N 
-ir  uai.  ot  *  -  I - — - - 

FORMAT  I'  FnTFR  IMF  OffiRFF  OF  FITTING  P 01  »N OM I A L * // » 

F  UkHAI «  > 

CALL  ORTHLSI T IMF  ?.P? .U .N? ,0 .0 »COI F .Al PHA .BEIA.K.ITT1. 
TTI4.4TT4.1N0U  - 

CAU  COF FSlO.rOFF.At PHA.PrTA.A.A.TIT 1 .1 1  T?  .1  T  T  T.  1  ND?  I 

. — -N4~=-K-  *  1- - -  -  - - —  — - 

UR  IT  F  (6  •  99  I 

GO  I  HOI  I  ^  141 

uRiTrih.inm  i.aiti 

CONI INUF 

LOROR  :  II 

NO  4000  - 1  —  1  .N4  - - 

PP  r  n 

A  ;  T  IMF?!  It 

no  ?nm  J  :  i  .a i 

PP  :  PP  •  AL4l****!4-ll 
CONTINUF 

ERROR  -  ERROR  *  CRRTI  »-RB  f»*R - 

CONT INUf 

LRROR  -  ERROR/N? 

ERROR  =  SORT  IFRRORI 
UR  ITF  (6.  TOO  I  roaoa 

FORMAT!’  STANOARO  OFXlATTON  OF  PRFSSURF  =  ‘.GPO.*//! 

- EOBAeUU.'P  .ISA.  *  A I II  *-441 - - - 

FORM  AT  ISA. I I.10X.GIS.7I 
GO  TO  RHP 
Rf  TURN 
F  NO 


- ,  !*4». 


TEST  RUN  :  XcP/CTO  CALIBRATION  PROGPAN  -HP  SAUnOFRS  -  IA733  -NOBOAJJJ 
iPRT.S  M.PLOTfR/XCP-CALIB 


TEST  RUN  :  XCP/CTO  CALIBRATION  PROG  PAN  -  K  It  SAUNOFRS  -  M7JJ  -NOROA33I 


COOC  33  IbUORKF  I  LEI 
I 

7 

3  - 

A 

5 

6 
7 
P 

A  -  - - 

10 
II 
17 
1  3 
I  A 

-15 - - 

16 
17 
If 
19 

?n 

71 

77 
73 

2A 

75 

76 

-  -7-7- 

78 

79 

30 

31 

37 
33 

3A 

35 

36 

37 
3  8 
39 
AO 
A  1 
A? 

A3 
AA 
AS 
A6 
A  7 
A  8 
A  9 

50 

51 
57 
S3 
5  A 

55 

56 


1  I  .PLOTER  /XCP  -CAL  I  B 

SUBROUTINE  PLOTER 
CANNON  /O! AGS /  IOIAGIIOI 


" C88AAA 
CONN  ON 
CONN  ON 
CONN  ON 
CANNON 


/ ROOC 1 /  lOOCRIIOOl 
/ ROOCF  /  F 00 CRT  TOO  I 
/ROOC A/  A00CRI100) 
/RHOB/  I R MORI  TOO 1 


CONNON  /BLOCK]/  TIN  EX  1 10001 .TX  I  1  000 1 »PI I  10001 .TC I  10  001 .PCI  10301 

COMMON  /BLOCK]/  TEMPI  ICOfll  .P«  10001.  TINE  I  10001 
CANNON  /PlfiCKA/TfNP?(  lOOOT  ^7310603  .TiPf  74  18001 
CONNON  /flLOCKS/T  1.T7.0T .NT .N? 

f>E AC  N<  IPflOl.PNT  IOCW».VVHUflf>».T?«10BP» 

f  OUI  VALENCE  INI  1  I.VI  11 1  .  IPNI1  I  .VRIlOnill  .(UK41I.VRI  70fl  1  I  I 


c 

c 

c 

c 

€-  -  - 

c 

c 

c 


inon 

c 

c 

c - 

c 

c 

c 

c 

c 

c 


SPFFOS  ANO  A  5  SBC  FA  ECU PRESSURES 


00 


limn  1  r  ?  .N7 

-U- N-7—  -4  -♦  -7- - 


Win  ;  «P?|J|-P7«J-M  1/ITINE7IJI-TTNC7I  J-l  II 
P  Ni 1 1  —  IP? Ul«P7IJ-H »/?. 

TTII1  =.  T  I  *1.7  I  Jl 
CANT  I  HUE 


PLOT  THE  PRESSURE  TINE  DERIVATIVE  AS  A  FUNCTION 
OF  PRESSURE 


Nl  =  N7  -1 
CALL  CONPRS 

CALL  TITtFl'nP/OT  VS  PRf  SSURF  t  •»  100  »  ’OP/DT  I  ON/SEC  I* • .1  On. 
’PRESSURE  EASES’  .180.  7.  .9.1  -  - 

call  GRAFI0..7..1D..n..inn..600. I 
CAM.  frame 
CALL  CROSS 

CALL  CURVEIWI7I.PMI71.NI.N1 
CALL  FNDPLI1I 


snor 


no  snno  \  -  i  ,N7 

MR  ITFI  70. SOU  I.bill.PVlII 

continue 

no  6onn  i  -  ?  ,h? 

ww  tit  -C.7s a  iw  i  i-ii*vu*ni*n.5*Miii 


26 


I 

I 


K 


TfST  RUN 


xcp/cro  cii  ibb*  rloh  pRor, »»« 


K  D  SiUNDFRS 


x*7u  -*o»DMn 


57 
SB 
59 
60-  — 
61 
6? 
t>1 
69 

65 

66 
67 
6  B 

69 

70 

71 
7P 
75 
79 

75 

76 

77 
76 
79 
BP 
81 
8? 

8  5 

*9 - 


6000  COMTINUt 

till  ftSNPtt?! 

CiLl  TITlf  «  *t)P/OT  VS  PRf55URr  -  F  U  T  T  Pf  0  »•  ,1 00  . 

— - 1 - «f>Pift7  tflOFSFf  U*t»00. - - - - - 

1  •P9F5SURF  (OB  It* .100. 7. .9. I 
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